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Placing a quantum dot close to a metal nanowire leads to drastic changes in its radiative decay 
behavior because of evanescent couplings to surface plasmons. We show how two non-Markovian 
effects, band-edge and retardation, could be observed in such a system. Combined with a quantum 
dot p-i-n junction, these effects could be readout via current-noise measurements. We also discuss 
how these effects can occur in similar systems with restricted geometries, like phononic cavities and 
photonic crystal waveguides. This work links two previously separate topics: surface-plasmons and 
current-noise measurements. 
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I. INTRODUCTION AND MOTIVATION 

When a photon strikes a metal surface, a surface 
plasmon-polariton (a surface electromagnetic wave that 
is coupled to plasma oscillations) can be excited. The 
concept of plasmonics^, in analogy to photonics, has 
arisen as a new and exciting field since surface plas- 
mons reveal strong analogies to light propagation in con- 
ventional dielectric components^ and provide a possible 
miniaturization of existing photonic circuits^. 

In a related context, a complete understanding of the 
dynamics of quantum systems interacting with their sur- 
roundings has become desirable, particularly with respect 
to applications for quantum information science. While 
the Markovian approximation is widely adopted to treat 
decoherence and relaxation problems, the non-Markovian 
dynamics of qubit (two-level) systems have come under 
increased scrutiny^. This is because a simple Markovian 
description is not adequate when the qubit is strongly 
coupled to its environment. In solid state systems, an 
exciton in a quantum dot (QD) can be viewed as such a 
two-level system. Recently single-qubit gate operations 
on QD excitons have been studied experimentally^. Fur- 
thermore, with advances in fabrication technologies, it is 
now possible to embed QDs inside a p-i-n structure^, such 
that electrons and holes can be injected separately from 
opposite sides. This allows one to examine the exciton 
dynamics in a QD via electrical currents^. 

Motivated by these recent developments in plasmon- 
ics and quantum information science, we show in this 
work how non-Markovian interactions between QD exci- 
tons and nanowire surface plasmons give rise to two in- 



teresting effects: band-edge and retardation. In a differ- 
ent system, the band-edge effect was originally predicted 
using the isotropic band-edge model^: the quadratic dis- 
persion relation, cok = u;c-\-A{k — kc)'^^ leads to a photonic 
density of state p{uj) at a band-edge ujc^ which behaves 
as l/\/uJ — uoc for uo >uOc- In a nano-wire, the band-edge 
effect stems from the non-linear behavior of the plasmon 
dispersion relation, in which there are similar quadratic 
local extremes at certain frequencies. The other effect 
we investigate here, retardation, is the multiple time de- 
lay of emission and absorption of plasmons between two 
QDs. With the incorporation of the system inside a p-i-n 
junction, we show that both effects can be readout via 
current-noise measurements. The possibility of observ- 
ing such phenomena in a QD spin qubit confined in a 
phononic cavity or a QD in a photonic crystal waveguide 
are also discussed. 
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II. BAND-EDGE EFFECT 



Consider now a semiconductor QD near a cylindrical 
metallic (we will consider silver here) nanowire with ra- 
dius a and longitudinal axis z as shown in Fig. 1. The QD 
and nanowire are assumed to be separated by a dielectric 
layer^. The n-th surface plasmon mode's components of 
the electromagnetic field at the surface can be obtained 
by solving Maxwell's equations in a cylindrical geometry 
(p and if denote the radial and azimuthal coordinates, re- 
spectively) with appropriate boundary conditions'^. The 
dispersion relations of the surface plasmons can be ob- 
tained by numerically solving the following transcenden- 
tal equation: 
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FIG. 1: (Color online) Schematic view of the system: a metal- 
Uc (e.g., silver here) nano-wire is embedded inside a GaN ma- 
trix and a (blue) QD (quantum dot) is placed on top of it. 
An evanescent electromagnetic wave couples the metallic wire 
and the QD. The exciton in the QD (presented by the two 
disks) can recombine, spontaneously emitting photons (green 
arrow) that produce surface plasmons on the wire (illustrated 
by the surface effect). 
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whose solutions are the dispersion relations uOn = ujn{kz). 
Here, / (0) stands for the component inside (outside) the 
wire. Also, Jn{Kip) are Hn\Kop) are the Bessel and 
Hankel functions, respectively. The dielectric function is 
assumed as 
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where eoo = 9-6 (for Ag), eoo = 5.3 (for GaN), uop 
is the plasma frequency, and r is the relaxation time 
due to ohmic metal loss^^. The magnetic permeabilities 
III and iio are unity everywhere since here we consider 
nonmagnetic materials. The reason to choose a silver 
nanowire here is that the plasmon energy TiijOp of bulk sil- 
ver is 3.76 eV with the corresponding saturation energy 
?iix}p/\/2 ^ 2.66eV in the dispersion relation. As we shall 
see below, variations of the dispersion relations in energy 
just match the exciton bandgap of wide-band-gap nitride 
semiconductor QDs. In related work. Gallium nitride is 
used as a matrix interface between a silver film and a in- 
dium gallium nitride quantum well"*^^. This is primarily 
because the refractive index of GaN reduces the surface 
plasmon energy to match that of the exciton energy. 




FIG. 2: (Golor online) (a) Dispersion relations Re[cOn] versus 
Kz of surface plasmons for the first few modes (n = 0, 1, 2, 3). 
The units for the vertical and horizontal lines are Q — uj/ujp 
and K = kzc/oop. (b) The enlarged plot of the dispersion 
relations of (a) in the regime of small kz- (c) Correspond- 
ing (Markovian) spontaneous emission (SE) rates into surface 
plasmons. As seen here, the SE rates are greatly enhanced at 
certain values of cjq- (d) Non-Markovian decay dynamics of 
QD excitons for 6 = 0.2/3 (black line), 0.4/3 (dotted line), and 
0.8(3 (dashed line). When ^ = 0, the red curve represents the 
result of the contribution from the n — 1 mode. 



The dispersion relations for various modes n are shown 
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in Fig. 2(a) with effective radius R = 0.1. The unit of the 
effective radius R (= cOpa/c) is roughly equal to 53.8 nm. 
The behavior of the n = mode is very similar to the 
two-dimensional case^^, i.e. Q gradually saturates with 
increasing wave vector kz. This is because the fields for 
the n = mode are independent of the azimuthal angle 
(p. However, the behavior for the n ^ modes are quite 
different. The first interesting point are the discontinu- 
ities around uj/c ^ kz. Further analysis shows that the 
solutions of uj are "almost real"^^ when kz > Re[uo\/c. 
Thus, the first Hankel function of order n, Hn\K^p)^ 
decays exponentially. This means that the surface plas- 
mons in this regime are confined to the surface {hound 
modes). For kz < Re[uj\/c^ however, the solutions of uo 
are complex^ as shown by the dashed lines in Fig. 2(b). 
H^\k^P) in this case is like a traveling wave with finite 
lifetime {non-hound modes). 



A. Spontaneous emission rates 



is the detuning of the radiation mode frequency u;n{kz) 
from the excitonic resonant frequency cjq, and gn,kz — 

do' En,kz is the atomic field coupling. Here, do and En.kz 
denote the transition dipole moment of the exciton and 
the electric field, respectively. The subindex "ex-sp" in 
-f^ex-sp refers to excitons (ex) and surface plasmons (sp). 

Assuming that initially there is an exciton in the 
dot with no plasmon excitation in the wire, the time- 
dependant wavefunction of the system then has the form 



\^lj{t)) = be{t) IT, 0)^^bn,kAt) li, e- 



(5) 



n,kz 



The state vector ||, 0) describes an exciton in the dot and 
no plasmons present, whereas ||, ln,/e^) describes the ex- 
citon recombination and a surface plasmon emitted into 
mode kz. With the time- dependent Schrodinger equa- 
tion, the solution of the coeflicient be{t) in 2:-space is 
straightforwardly given by 



Once the electromagnetic fields are determined, the 
spontaneous emission (SE) rate, Tse, of the QD exci- 
tons into bound surface plasmons can be obtained via 
Fermi's golden rule. The SE rates of the first few modes 
(n = 0, 1, 2, 3) are shown in Fig. 2(c) with effective radii 
R = 0.1. In plotting the figures, the distance between 
the dot and the wire surface is fixed as d = 10.76 nm. 
The novel feature is that the SE rate approaches infin- 
ity at certain values of the exciton bandgap ujq. Mathe- 
matically, one might think that at these values the corre- 
sponding slopes of the dispersion relation are zero. Phys- 
ically, however, this infinite rate is not reasonable since 
it is based on perturbation theory. Therefore, one has 
to treat the dynamics of the exciton around these val- 
ues more carefully, i.e. the Markovian SE rate is not 
enough. One has to consider the non-Markovian be- 
havior around the band-edge, which means the band 
abruptly appears/disappears across certain values of oj. 



B. Non-Markovian dynamics 

To obtain the non-Markovian dynamics of the exciton, 
we first write down the Hamiltonian of the system in the 
interaction picture (within the rotating wave approxima- 
tion). 
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where aij = |i) (j|(i,j =T, i) are the atomic operators; 
an,kz ^ are the radiation field (surface plasmon) 

annihilation and creation operators; 

An,/c, = (^n{kz) -^0 (4) 
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Z + i[Un{kz) - UJo] 



(6) 

In principle, be{t) can be obtained by performing a nu- 
merical inverse Laplace transformation to Eq. (6). 

To grasp the main physics and without loss of gener- 
ality, we focus on the values of coq close to one particular 
local extremum (maximum or minimum), where the dis- 
persion relation for a particular mode becomes quadratic. 
In this case, the dispersion relation for this particular 
mode around the extreme can be approximated as 
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where the extremum is located at {kn^c^^n.c)- The +/— 
sign represents the approximate curve for the local mini- 
mum/ maximum of the dispersion relation. Once we make 
such an approximation, the radiative dynamics of the QD 
exciton is just like that of a two- level atom in a photonic 
crystal^' with 
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where 
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for local minima (8) 



for local maxima (9) 



(10) 



is the detuning to a specific extremum and 7 is the de- 
cay rate contributed from other modes. For example, 
hereafter we choose coo to be close to the minimum of 
the n = 1 mode, and thus only this n = 1, and the 
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n = mode, strongly interact with the exciton. The 
other modes can be treated as a (Markovian) decay pro- 
cess with a rate 7. 

The coefficient be{t) can now be obtained^' by per- 
forming the Laplace transformation to Eqs. (8,9). The 
black, dotted, and dashed lines in Fig. 2(d) represent 
the decay dynamics of the QD excitons for different de- 
tunings: S = 0.2/3, 0.4/3, 0.8/3, respectively. Here, /3 is 
the decay rate of the QD exciton in free space. As men- 
tioned above, when plotting Fig. 2(d), uoq was chosen to 
be close to the local minimum of the dispersion relation 
of the n = 1 mode. The radius of the wire and the wire- 
dot separation are identical to those in Fig. 2(a). As can 
be seen in Fig. 2(d), there exists oscillatory behavior in 
the decay proffie of \be{t) \ , demonstrating that the decay 
dynamics around the local extrema is non-Markovian. If 
one only considers the contribution from the n = 1 mode 
and set the detuning 5 = 0, the probability amplitude 
would saturate to a steady limit, as show by the top 
red curve in Fig. 2(d). This quasi- dressed state is rem- 
iniscent of damped Rabi oscillations in cavity quantum 
electrodynamics, and also appears in systems of photonic 
crystals^'^^. 
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FIG. 3: (Color online) A schematic diagram of a p-i-n junction 
with a quantum dot (QD) evanescent ly coupled to surface 
plasmons in a nanowire. 



measurements; i.e. by placing a QD p-i-n junction close 
to the nanowire as shown in Fig. 3. 

In addition to the Hamiltonian i^ex-sp in Eq. (4), we 
now need to consider the tunnel-couplings to the electron 
and hole reservoirs^: 



E 

q 



(v^io)(Ti + w^q4io)ai+if.c.), (11) 



where Cq and are the electron operators in the right 
and left reservoirs, respectively. Here, Fq and Wq couple 
the channels q of the electron and the hole reservoirs. 
We also introduced the three dot states: |0) = |0,/i). 
It) = |e,/i), and ||) = |0,0), where |0, /i) means that 
there is one hole in the QD, |e,/i) is the exciton state, 
and |0,0) represents the ground state with no hole and 
no excited electron in the QD^. 

Together with Eq. (3), one can now write down the 
equation of motion for the reduced density operator 



Jo 

[HT{t') + H,^_,p{t')M')]], 



(12) 



where S(t') is the total density operator. Note that 
the trace, Tr, in Eq. (12) is taken with respect to both 
plasmon and electronic reservoirs. Without making the 
Markovian approximation to the exciton-plasmon cou- 
plings, one can derive the equations of motions of the 
dot operators^^. With the help of counting statistics, the 
noise spectrum is then given by 



where 
Bico) = 



Si^{u;)=2eI{l + rR[B{io) + B{-io)]}, (13) 
A{iuj)rL 



A{iLo)rLTR + {A{iLo) + iLo){rL + iLo)(TR + ito) ' 

(14) 

Here, / is the stationary current, Tl and F^^ are the 
tunneling rates from the electron and hole reservoirs, and 
A{z) = c{z) + c*(2;), where 



c{z) = E 



Z + i\uJn(kz) - ^o] 



(15) 



Readout of the band-edge effect via 
current-noise 



With recent advances in fabrication technologies, it is 
now possible to embed QDs inside a p-i-n structure^. 
Furthermore, the interest in measurements of shot-noise 
in quantum transport has grown recently owing to the 
possibility of extracting valuable information not avail- 
able in conventional dc transport experiments^^. We thus 
propose to bring these two branches of condensed mat- 
ter physics together: surface-plasmon and current-noise 



Figure 4(a) shows the noise spectrum Si^{uj) as a func- 
tion oibj. As for Fig. 2(d), the value of cjq here is chosen 
to be close to the local minimum of the n — \ mode. 
Here, the F^ and F^^ are set equal to 0.01/3 and 0.1/3, re- 
spectively. The solid (dashed) line represents the result 
for (5 = cjo — ^n,c = —0.01/3 (0.01/3). The interesting fea- 
ture here is that there are discontinuities at = ±0.01/3. 
For the case of b = —0.01/3, the Poissonian value of the 
noise spectrum \Si^{iS) = 1 for -0.01/3 < uj < 0.01/3] 
is analogous to that of putting a two-level emitter inside 
the band-gap, while, for S = 0.01/3, the sub-Poissonian 
value is the situation outside the band-gap^^. Figure 4(b) 
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FIG. 4: (Color online) (a) Noise spectrum as a function of ou. 
Like that in Fig. 2(d), the value of ooo here is chosen to be 
close to the local minimum of the n — 1 mode. Here, Fl and 
Fr are set equal to 0.01/3 and 0.1/3, respectively. The solid 
(dashed) line represents the result for S = —0.01/3 (0.01/3). 
(b) Density plot of the current-noise spectrum as functions of 
00 and detuning 6 = ooo — C0n,c, which are both in units of /3, 
the decay rate of the QD in free space. As seen here, there 
are discontinuities along the lines 6 — zbcj, which is an indi- 
cation of the band-edge effect. In the inset we show the same 
calculation using the numerical solution for the dispersion re- 
lation, not the quadratic approximation. This illustrates that 
the important features are all contained within the quadratic 
approximation. 



shows the density plot of the noise spectrum as functions 
of both u and detuning, S = ujq — ^n,cAs seen there, for 
^ < 0, the values of Si^{u;) in the regime — \S\ < uo < \5\ 
are larger than those in a; < and \5\ < co. For 

J > 0, however, it is the opposite behavior. In addition, 
one also observes that there are discontinuities along the 
lines S = Together with the results in Fig. 4(a), 

we conclude that the feature of discontinuities in these 
noise spectra can actually be viewed as an indication of 



a band-edge effect. 



III. RETARDATION EFFECT 

By placing two QDs close to the nanowire, and by mak- 
ing use of the one-dimensional propagating feature of the 
nanowire surface plasmons, another non-Mar kovian ef- 
fect, the retardation, can be observed. For simplicity, 
the exciton energy hcoo of the two identical dots is set 
well below the local minimum of the n = 1 mode, such 
that only the n = mode contributes to the decay rate. 
Thus, the interaction Hamiltonian can be expressed as 



Hi 



-ih 



E 



kz ^ 



ikz zi 



IT)nUI).(16) 



where zi is the position of the l-ih. dot, and the distance 
of the two dots to the wire surface is the same. Assuming 
that only dot-1 is initially excited, the state vector of the 
system can be written as 

\m) = hit) \n,o)+b2{t) UT,o)+^6fe.(t) lu.in.fej , 

(17) 

with the initial conditions: bi{0) = 1, 62(0) = 0, and 
bkzi^) — 0- The time-dependent solutions are straight- 
forwardly given by ^1(2) (^) = [C-\-{t) ±C-{t)]/2 with 



-| pioo-\-e 

c±{t) = — ds- 

2in i_oo+e t 



(18) 



where 



G{s) = {5 + i[uJn(k^) - ojq]}-^ + {5 + i[ujn(k^) + uJo\}~\l^) 

Following the well-known treatment of retardation^^, one 
can obtain the probability amplitudes of the dots in the 
regime of /cqt > 3 



^1(2) W 



CO 



m=0,2,4 ... 
(m=l,3,5 ...) 
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ml 
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X 



where r = \z2 — z\\^ ko ^ loq/v^ v is the velocity of the 
surface plasmon on the wire, 70 is the spontaneous emis- 
sion (SE) rate of a single QD exciton into a surface plas- 
mon, and H is the unit step function. 

One might argue that the surface plasmons inevitably 
experience losses as they propagate along the nanowire, 
which could limit the feasibility of observing the retarda- 
tion effect. One solution to this would be to couple two 
QDs to two separate nanowires. Meanwhile, the wires 
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Si^iu) = 2eI{l^TR[nR{s = -iuj)^nR{s = iu)]}. (21) 

In Eq. (21), nR{s) is the Laplace-transformation of the 
ground state occupation probabihty ni^(t) = (|||) (|||)^, 
where the average is over both the electronic and pho- 
tonic reservoirs. The red and blue lines in Fig. 5(b) 
represent the noise spectra for 7or/'y = 27r and 47r, re- 
spectively. As seen there, the main difference to the 
non-retarded case (black line) is the oscillatory behav- 
ior, which depends on the inter-dot separation r. One 
recalls that in the non-retarded situation there should 
be no difference whenever ^qv/v = myr, where m is an 
integer^^. The green line in the inset of Fig. 5(b) is the 
result for ^qv/v = 1.9. This means that even if the value 
of 7or/'U is not equal to rnvr, one still could observe the 
predicted oscillatory behavior. 



IV. BAND EDGE EFFECT IN PHONON 
CAVITIES 

The non-Markovian effects studied above can also be 
observed in other physical systems. For example, if one 
considers a free standing slab^^ with width small elas- 
tic vibrations of a solid slab can then be defined by a vec- 
tor of relative displacement u(r,t). Under the isotropic 
elastic continuum approximation, the displacement field 
u obeys the equation 



= c?V^u+(c?-c?)V(V-u), 



(22) 



FIG. 5: (Color online) (a) Proposed device for the observa- 
tion of retardation effects via current-noise. The two QDs 
are coupled to two separate nanowires. Meanwhile, the wires 
are evanescently coupled to a phase- matched dielectric waveg- 
uide, (b) Current noises of the double-dot device. The red 
and blue lines represent the results for r/v = 27v and 47r, 
respectively. Recall that r is the inter-dot separation. The 
black line is the result for the Markovian case. The most ob- 
vious feature for the non-Markovian effect is the oscillatory 
behavior (red and blue curves). Inset: Noise spectra with 
(green line) and without (black line) retardation effects when 
r/v 1.9. 



where ci and Ct are the velocities of longitudinal and 
transverse bulk acoustic waves. To define a system of 
confined modes, Eq. (22) is complemented by the bound- 
ary conditions at the slab surface z = ±w/2. Because of 
the confinement, phonons will be quantized in subbands. 
For each in-plane component qy of the in-plane wave vec- 
tor there are infinitely many subbands. Since two types 
of velocities of sound exist in the elastic medium, there 
are also two transverse wavevectors qi and qt. If one fur- 
ther considers the deformation potential only, then there 
are two main confined acoustic modes: dilatational waves 
and flexural waves. For dilatational waves, the parame- 
ters qi^n and qt^n can be determined from the Rayleigh- 
Lamb equation 



would be evanescently coupled to a phase-matched di- 
electric waveguide^^. In this case, one could have both 
the advantages of strong coupling from the surface plas- 
mons and also long-distance transport in the dielectric 
waveguide. In addition, the non-Markovian retarding ef- 
fect can also be measured via current noise if one of the 
dots is embedded inside a p-i-n junction as shown in Fig. 
5(a). Following the procedure described above, the noise 
spectrum is given by 



tan (gt,n^/2) _ 
tan {qi^riw/2) 

with the dispersion relation 



M\\Ql,nqt,n 



(23) 



(24) 



where 



is the frequency of the dilatational wave in 



mode (n,q||). For the antisymmetric flexural waves, the 
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solutions qi^n and gt,n can also be determined by solving 
the equation 

tan _ _ M\\qi,nqt,n , . 

tanfe,W2)" {qI-qW ^ ^ 

together with the dispersion relation, Eq. (24). 

Figures 6(a) and (b) numerically show the dispersion 
relations for dilatational and fiexural waves, respectively. 
As can be seen in the insets, local minima also appear in 
the dispersion relations. An enhanced relaxation rate due 
to the phonon van Hove singularities has been predicted 
if a double-dot charge qubit^^ or single-dot spin state^^ 
is embedded in such a phonon cavity. However, as we 
have mentioned above, the greatly enhanced rates are 
also from the band-edge like effect^^, and one should treat 
the dynamics of the qubits as non-Markovian. As for the 
retardation effect, the two QDs may also be embedded 
inside a well-designed photonic crystal waveguide^^, in 
which the propagation of the photon is restricted to one 
dimension. In this case, the advantage of the retardation 
effect in one dimension is still kept, and the combination 
with the p-i-n junction should also be workable^^. 

V. CONCLUSIONS 

In summary, we have shown that nanowire surface 
plasmons, which we consider as a bosonic reservoir with a 
restricted geometry, have a non-linear dispersion relation 
with extreme values at certain frequencies. When cou- 
pled to a QD exciton (combined with a p-i-n junction) we 
described how it should be possible to observe the non- 
Markovian dynamics of these effects when the recombi- 
nation energy of the exciton is close to the bandgap of the 
plasmon reservoir. We calculated specific results for the 
current-noise frequency spectrum and observed unique 
signatures of these 'band-edge' non-Markovian dynam- 
ics. 

Furthermore, we have shown that the retardation ef- 
fect, another non-Markovian effect which occurs when 
two dots are both strongly coupled to the same nano- 
wire, has also unique signatures in the current-noise. Fi- 
nally, we illustrated how these effects might also be ob- 
served in a QD spin qubit (or double-dot charge qubit) 
embedded inside a phonon cavity. 
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